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Abstract 

The application of a modified gluon propagator and an infrared frozen running coupling constant 
in the description of the proton-proton elastic scattering is investigated. The model based on the 
lO ' exchange of the Pomeron as a pair of non-perturbative gluons is modified to include the frozen 

o : 

. coupling constant in the infrared kinematic region. A good fit to the data of the differential cross 

5^ \ section for different center of mass energies (53 GeV and 1.8 TeV) is obtained. 
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I. INTRODUCTION 



The Pomeron has an important role in the elastic scattering, confirmed by the description 
of the experimental data of the hadronic collisions and diffractive dissociation [l, Q, O]- In 
the framework of the Regge theory [3] , the Pomeron is the exchange of the vacuum quantum 
numbers, with the scattering amplitude given by 

A{s, t) oc = a(0) + a'{t) (1) 

The most simple picture within Quantum Chromodynamics (QCD) for the Pomeron is 
the exchange of a pair of gluons, as showed by Low 4] and Nussinov 5]. However, the 
applicability of perturbative QCD (pQCD) in the description of soft processes is not free of 
troubles, specially in the singularity structure of the scattering amplitude in the kinematic 
regime of low momentum transfer. The singularities in this regime include a large value of 
the strong coupling constant and collinear divergences in the gluon and quark propagators. 
Usually, when the standard gluon propagator is employed, a divergence (pole) exists when 
the momentum transfer goes to zero. A solution to this difficulty is the introduction of a 
cut by hand or the use of the spontaneous symmetry breaking mechanism to generate a 
mass to the gluon, with the price of the introduction of the Higgs bosons to the theory. 
Although a bare mass for the gluon (corresponding to a massive propagator) has problems 
with unitarity and renormalization, this propagator has some interesting phenomenological 
features l6|. 

Landshoff and Nachtmann (LN) [7] proposed a phenomenological model to overcome these 
difficulties, arguing that in this kinematic region the non-perturbative effects are important. 
In this model, these effects change the gluon propagator, which is related with the QCD 
vacuum through the gluon condensate 

Ross use the ideas of the LN model in a phenomenological model, in which the Pomeron 
exchange is described as in the Balitski-Fadin-Kuraev-Lipatov (BFKL) approach, mean- 
ing a gluon ladder with the exception of the gluon propagator, modified in the infrared 
region. The correct model must include the infrared modification of all the Green functions 
as quark propagator and vertex parts, but the hope is that all infrared contributions can be 
included in the gluon propagator. These ideas are employed by Cudell and Ross in the 
description of the pp elastic scattering in low momentum transfer. 
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In this paper, we propose a modification in the above modeL Despite the fact that a 
running couphng is a next leading order contribution, we will use a dynamical coupling 
with a frozen (finite) infrared behavior, but related with the gluon propagator as analyzed 
in previous works [3, U, ll^ to describe the process of pp elastic scattering at different 
energies. 

The paper is organized as follows. In the section m the previous models of the pp elastic 
scattering are discussed. In section IIIIl a brief review of the results of gluon propagators 
and frozen running couplings are presented. The modified phenomenological model and 
its results are presented and discussed in the section HYl The conclusions are presented in 
section El 

II. THE MODEL 

The model employed to describe the proton-proton elastic scattering in this work was 
proposed by Cudell and Ross , based on the two-gluon model to the Pomeron. A remark- 
able feature of this model is the strong dependence on the wave functions of the involved 
hadrons, giving a possible source of uncertainty. In the above work, the LN model Q for the 
Pomeron is used to avoid the infrared divergences employing a gluon propagator calculated 
from Dyson- Schwinger equations as shown in the next section. 

The hadron-hadron scattering can be derived from the photon-photon scattering with 
appropriated form factors. With the hadron constituted by valence quarks, the scattering 
amplitude for two gluon exchange can be written as [l^ 

A^^'''' = ^sal n^n^ C j dKdK 6^^^ (A - k„ - k,) V{kl)V{kl) x 

[^i"^(k, + k,) -4^(k„k,)] [S'^'{K + kb)-S^%kaM)] , (2) 

where C is the color factor of the process (C = 8/9), is the number of valence quarks 
{nafi = 3) and hi is the type of hadron {ha^h = p). After carrying out the integral, the 
resulting expression is 

Af = ^isal jdkV(k'')V{{A-kf)[Si{A)-S2{k,A-k)]\ (3) 

The form factor Si can be identified with the elastic Dirac form factor for the proton, 

4m2-2.79t 

Si{t = -A^ = ^ (4 

^ ' (4m2-t)(l-t/0.71)" 



where rrip is the proton mass. The second form factor has a more arbitrary structure, being 
fixed requiring that infrared divergences disappear. A simple form that agrees with this 
assumption is 

£2{K, k,) = Si{kl + kl- fK ■ kfe), (5) 

where the / parameter varyes depending on the specific hadron. In the proton case, its value 
is choosen such that the quark wave functions are peeked at /3 = 1/3. In ^], / = 7, but 
using other arguments [15] the value of / = 1 is also advocated. 
The scattering amplitude can be rewritten as 



APP{s,t) = Sisal 



where 



r, = |d^kp(| + k)l)(|-k)G^(g,0), 



(6) 



(7a) 



q 



q 



d'k P + kj I? - kj Gp [q, k - I j [2G,(g, 0) -Gp^q,k-^)\, (7b) 

where s,t = —q^ are the Mandelstam variables and Gp(g, k) is a convolution of the proton 
wave functions, 

G,(,.*=) = /d^pda#(a.p)V.(a.p-*-a,). 
and can be related with the Dirac form factor as follows, 



Gp(g,0) = Fi(g2) 



k' 
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(8a) 
(8b) 



The terms in eqs. fl7a|7bj) have a simple interpretation: the former comes from Feynman 
diagrams where the gluons are connected to the same quark in the proton. The last one 
comes from diagrams where the gluons are connected to different quarks in the proton j^. 

The total cross section is given by the optical theorem, a^^^ = A2^{s,0)/is, as well as 
the differential cross section da^/dt = |^2^(s, t)p/167rs^. For low momentum transfer, the 
elastic differential cross section is fitted by a exponential expression, da/dt = Ae^^ = 
^tot e'^V(16vr), where B is the logarithmic slope, given by i? = d/dt\ln{da/dt)]\t=o- 

With the scattering amplitude above, a finite result for the total cross section is obtained 
due to the cancelation of the infrared divergences 16|. Otherwise, the differential cross 
section has an unphysical result {B{t = 0) ^ oo) if the perturbative gluon propagator is 



employed. When the Landshoff-Nachtmann model is used, the contribution from the term 
7^ is negligible and the remaining infrared divergences are regularized by the modified gluon 
propagators. 

However, the two gluon description of the scattering amplitude is a crude approximation, 
which does not yield any energy dependence in the above cross sections. The reason is the 
absence of the gluon ladders in the s channel, as in the BFKL framework, which gives an 
energy dependence in the cross sections. To restore this dependence, we assume that the 
amplitude above only describes the energy independent part of the cross section and an 
extra term, introduced ad hoc gives the Regge energy behavior, namely, 

where ap(t) = a{0) + a'{t) is the (soft) Pomeron trajectory (with a{0) = 1.08 and a'(t) ~ 
0.25 GeV~^ [17]) and Sq is an energy scale. The extra Regge term gives the following cross 
sections, 



ap(0)-l 

So, 



c^tot = ( — ) o^tot, (9a) 



da f sV^-^'^-' da^ 

(9b) 



dt \Sq J dt 

The above model was used by 0, 16, Q where the non-perturbative propagator comes 
from different methods, discussed in the next session. The experimental data are ISR results 
for elastic proton-proton scattering at = 53 GeV jl9l |. 



III. GLUON PROPAGATOR AND RUNNING COUPLING CONSTANT 

A central idea of the LN model for the Pomeron is the infrared modified gluon propagator, 
where the modifications are induced by the QCD vacuum effects, and its non perturbative 
character required the application of methods of the same nature to obtain the gluon prop- 
agator. The most popular methods are the Dyson- Schwinger equations (DSE) [20|, |21[ and 



the numerical computational simulations on lattice field theory 

The DSE's are an infinite system of non-linear, coupled integral equations which relate 
the different Green functions of a quantum field theory. In the case of QCD, due to the 
complexity of the system to be solved, there are several analytical approximations in the 
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literature to obtain a solution. A discussion among the different methods and approximations 



employed can be found in 



20 



2l|. 



The lattice field theory, in opposition to the DSE, is a powerful pure numerical method, 
where the main idea is that the divergences of the theory are regularized by the discretization 
of the space-time. Otherwise, the lattice approximation is not free of problems: finite size 
of the lattice, spacing between the sites of the lattice and fermion simulation. 

Due to the approximations employed in the resolution of the problems pointed above, 
there are many different solutions for both methods in the literature. We focus here on the 
solutions which have the particular property of a dynamical gluon mass. A gluon propagator 
with this feature was used in successful descriptions of several processes |lO, luJ, • The 
first solution for the gluon propagator which will be used in this work is the Cornwall's 
solution [23I, from DSE's in the axial gauge with the use of a technique of resummation of 
Feynman diagrams. 



with 



r 2 I 2/2 



hg^ In 



q^ + 4m^(q' 



^^QCD 





-12/11 



(10) 



(11) 



g - 500 ± 200 MeV for Aqcd = 300 MeV and h = 33/(487r2) is the leading order 



where m 

coefficient of the (3 function. A variation of the above propagator is given by 



Other solution is obtained by Habel et al. 
which gives the following propagator. 



(12) 



l25l | using another type of approximation. 



q24-64' 

where 6 is a free parameter. Gorbar and Natale [2q use the propagator 



"' = q^ + /i?e(eV?-q^) + ^e(q^ 



(13) 



(14) 



to calculate the vacuum QCD energy through effective potentials, where 6 is the step 
function, Hg is related with the gluon condensate and ^' is a calculated parameter. An 
alternative form used in this work is 



[^Giv(q')]"' = q' + "^GAre(m: 



GN 



q^) + 



m 



GN 



q- 



e(q^ 



(15) 



where is a mass parameter. A recent result for the gluon propagator, obtained from 

the DSE's in the Mandelstam approximation, is due to Aguilar and Natale [27 1, where the 
gluon propagator is well fitted by the expression 



q2 + A^2(q2)' 



rrin 



(16) 



Another recent result is the solution of DSE's in the Landau gauge calculated by Alkofer 
1 2^, 1^ , where the gluon propagator are given by 



and collaborations 
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UJ 
~2 



q 



q- 



^^■QCD + q 



2k 



(17) 



where uj = 2.5, Aqcd 
coupling constant, 

1 



510 MeV, K ^ 0.595, 7 = -13/22 and ai^^\ci^) is the running 



(qVA2 



QCD) 



«s(0) + 



4:71 q2 



+ 



(qVA' 



QCD 7 



(18) 



where a,(0) ^ 2.972 and po = 11. 

The momentum dependences of the gluon propagators are displayed in the figure d 
in comparison with the perturbative one. The most remarkable feature is the infrared 
behavior (k ^ 0), where exists a finite value in the case of the Cornwall propagator, for 
example; a zero value, as in the Habel/Gribov propagator, and the divergent perturbative 
propagator. Another remarkable fact is the same ultraviolet behavior (k 00) of the 
modified propagators and the perturbative one. 

Nevertheless, there is a narrow relation between the gluon propagator and the running 
coupling, widely discussed in [l^ [ll | . The phenomeno logical applications of a frozen cou- 



10, 



12 1 and other 



pling constant include heavy quarkonia decays j31|, meson form factors 
observables |llj]. In some of the above results (specially in Q), the data available are de- 
scribed when the modified gluon propagator and the frozen coupling constant are used at 
the same time. 

The analytical expression for the running coupling is obtained following the Cornwall's 
solution j23,|3^, namely. 



a,(A;2) =47r / PqIu 



^^QCD 



(19) 
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FIG. 1: (a) Distinct modified gluon propagators compared with the perturbative one. The propa- 
gators are indicated in the plot legend, (b) Comparison between different results for the running 
coupling constant. In the detail, two results of the Cornwall frozen coupling with two different ^'s. 
The parameters of the curves are given in the text. 

where the explicit momentum dependence is the same as the one loop perturbative calcula- 
tion, but with an extra term, ^m^(/c^), resulting in the frozen infrared behavior where ^ is 
a parameter that varies with the propagator and measures the frozeness of the coupling. In 



the original work j2c 



22], the value of this parameter is obtained as ^ ~ 4 or larger. The 
massive term, mp(fc^), also depends on the gluon propagator employed, being the dynamical 
mass term of the propagator. For example, in the case of the Aguilar and Natale propagator, 
eq. this term is m-p{k'^) = Ai'^{k'^). The exception is the Alkofer case, with a proper 

expression for the running coupling constant, displayed above. 

In the fig. ^ we present the different a^'s (given by eq.()19|)). where the freezing of the 
running coupling constant, the matching of the frozen running coupling, and the behavior 
of the perturbative result in high momentum are very clear. The last result comes from the 
behavior of the massive term (m-p(A;^)) which goes to zero when the momentum is large in all 
the propagators employed in this work (in the case of massive propagator, the mass is small). 
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The parameters of some propagators (m^, b, mo,. . . ) are determined when the results of the 
elastic pp scattering are fitted (see below) with exception of the Alkofer propagator. The 
figure also shows the ^ parameter dependence of the Cornwall solution, eq. (fTIH) . As seeing 
in fig. n and also in the results of the model, there is a small difference for distinct ^ values. 



IV. THE MODIFIED MODEL 



In the previous model for the scattering amplitude, the running coupling constant was 
consider fixed. In the value of the coupling was determined throught the effective 
Pomeron coupling with the hadrons, which depends on the gluon propagator and has well 
defined experimental value. 

In this work, motivated by the good results of the previous calculations with a gluon 
propagator and a frozen running coupling |12], we modify the previous model of pp scatter- 
ing, considering now a frozen running coupling constant and a modified gluon propagator. 
We emphasize that this is a phenomenological modification, since a running coupling is a 
next leading order contribution and the model employed is a first order contribution to 
the Pomeron exchange. As usual in the previous works, the scale of the running coupling 
constant is chosen as the incoming momentum into the quark-gluon vertex. 

The scattering amplitude in the modified model then reads 

gM2 



G,{q,0)-G,[q,k-^j\ . (20) 

where the terms from eqs. (fTaj) and (|7bj) were combined as one term only and the running 
coupling terms are included in the integrand. 

The differential cross section is then calculated using the eq. (pUj) into eq. (|9E|). varying 
the different propagators as well as the running couplings constants, as pointed out above. 
The results of the calculation are displayed in the fig. El for two sets of experimental data 
from different center of mass energies: y/s = 53 GeV from ISR and y/s = 1.8 TeV from 
Tevatron [s^ . 

We obtain a good agreement with the data in both energies and in the displayed range 
of transfered momentum. The slope of the curves are in accordance with data, due the ad 
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FIG. 2: Results of modified model, fitted to the pp elastic scattering data with ^/s = 53 
GeV (right) and ^ = 1.8 TeV 3( 



I (left) with different gluon propagator and frozen coupling 



constants. The parameters of the fits are given in the text. 

hoc term in the amplitude. The energy independent part of the amplitude, modeled by the 
eq. © gives a contribution on the t behavior and the intercept {da / dt\t=o) of the curve. 

Some considerations about the results can be made. First, the variation of the Cornwall 
propagator, eq. (fT^ . and the Gorbar/Natale one, eq. (fT^ are not displayed in the fig. El 
since these propagators give results quite similar to the presented ones. 

The values of the gluon propagator parameters used to obtain the results shown (or not) 
in the figure are: Cornwall, rrig = 0.53 GeV (with Aqcd = 0.3); Cornwall variant, rrig = 0.3 
GeV; Gorbar/Natale, rrig = 0.4; Habel/Gribov, b = 0.3 and Aguilar/Natale, rriQ = 0.35. The 
calculation is sensitive to the value of the parameters, giving larger modifications when the 
massive parameter is changed. In opposition, when the ^ parameter is changed, the change 
in the result is significantly smaller than with the mass parameter. The reason is the small 
modification in the running coupling when ^ value is changed, as pointed out in the previous 
section. Thus, we maintain the value of the ^ fixed in 4. The results of the modified model 
show that the mass parameter, for this process, is in the range of 300 MeV to 600 MeV, in 
agreement with the previous works which employ the modified gluon propagators. 
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When the Alkofer et al. propagator, eq. (fTTj) . and the corresponding running couphng 
constant, eq. (fTHj) . with the original set of parameters is employed, we obtain a result that 
overestimates the experimental data in both energies. Similar result was found by |lfll . 11 1 
for another observables. 



V. CONCLUSIONS 

In this work, we consider a frozen running coupling instead of a fixed one in a model 
for the amplitude of pp elastic scattering with the exchange of two non-perturbative gluons. 
The same idea was employed in another processes^which a small momentum transfer as, for 



lOL ll^ . with successful description of the 



example, the description of meson form factors 
data. 

In the present process, we obtain a quite good agreement with experimental data for 
distinct energies. The results obtained are along in the same lines of the previous attempts 
to describe the process, but considering other propagators. The results show that the employ 
of modified gluon propagators is valuable in phenomenology, but with some caution, since 
the model used to construct the scattering amplitude is a first approximation. A possible 
extension is to employ a complete set of Green functions: the quark and gluon propagator 
and the quark-gluon vertex. 

The application of the distinct results for the gluon propagator (and Green functions) 
found in the literature is also a field to test these solutions as well as the frozen running 
coupling constant. In the present work, the results point out to a gluon propagator with 



dynamical mass and frozen running coup 



dition with the previous ones 
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12 



ing constant, although the present results in ad- 
la], do not determine which is the most suitable 



expression for the propagator, mainly due the approximations employed and the still low 
knowledge in the literature on the interplay of whole aspects involved in these calculations. 
The consistency of these ideas used here can be checked in another type of diffractive pro- 
cesses, for example, diffractive photo-production of light vector mesons, which should deserve 
some analysis in this framework. 
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